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^ I We study large scale collective dynamics of isoscalar type and examine the in- 

fluence of interactions residual to independent particle motion. It is argued that 
for excitations which commonly are present in experimental situations such inter- 
actions must not be neglected. They even help to justify better the assumption 
of locality, both in time as well as in phase space, which is necessary not only for 
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^ ■ such classic approaches to collective motion as the "cranking model" but also for 

the more general picture of a transport process. With respect to dissipation, our 
results are contrasted with those of wall friction. 
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1 Introduction 

After the discovery of the shell model it has become customary to base the descrip- 
tion of collective motion on the picture of single particles moving independently within a 
deformed mean field. This approach was introduced in the early 50'ties by A. Bohr and 
B. Mottelson to portray low energetic collective excitations, and to the present day there 
can be little doubt that this approximation is adequate for that regime. It is somewhat as- 
tonishing, however, that this picture still is vindicated by many groups even for situations 
where the nucleons are heated up to considerable amount, say to temperatures of a few 
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MeV. After all, in the very early days of nuclear physics collective motion of large scale 
was considered to be governed by dissipative processes, which in turn imply the presence 
of fluctuating forces. Such a picture may be condensed into the one equation, which was 
suggested by Kramers [1] already in 1940 to describe nuclear fission. It reads 



d P d fdV(Q)\ d P ^ 



dQM dP \ dQ J dPM ' dPdP 



f{Q,P,t) (1) 



and has the structure typical of a Fokker-Planck equation. Actually, Kramers consider- 
ations strictly adhered to classical motion, in which case the diffusion coefficient is given 
by the Einstein relation Dpp — 7T; in the quantum case another term appears. As is 
well known, Kramers has used this equation to calculate the decay rate for a meta-stable 
situation like fission, in generalization of the famous Bohr- Wheeler formula. In these days 
the origin of dissipation was attributed to the strong "correlations" among the nucleons, 
as they can be understood within or follow from N. Bohr's compound nucleus — and which 
by definition occur "incoherently". 

In this lecture we want to look at this transition from " independent particle motion to 
coUisional dominance" in the view of the " linear response approach" , a complete version of 
which can be found in [2] . This discussion will be complemented by presenting new aspects 
in the relation to wall friction, following the more recent considerations of the group of 
W.J. Swiatecki, J. Blocki and others (see [3]). The applicability of linear response theory 
may be understood by the following arguments. First, one my note that the solution of 
(1) can be written in the following way 

f{Q, P,t) = j dgodPo K{Q, P, t; Qo, Po, to) f{Qo, Po, to) (2) 

where K{Q, P, t; Qq, Pq, to) is interpreted as the conditional probability for the system to 
move from Qq, Pq at to to Q, P at time t. On both sides of this relation the distribution /, 
the "joint probability" , may be replaced by conditional probabilities defining the transition 
say from a to to the final time t through an intermediate step at ti. The resulting relation 
is nothing else but the Chapman-Kolmogorov equation. (For a discussion of such general 
properties we may refer to the book by van Kampen [4]). The procedure just described 
may be repeated as often as one likes. Starting from the given equation (1) one may 
introduce arbitrarily small time steps St in completely rigorous manner. The reason for 
this behavior is found in the fact that this equation (assumedly) describes a genuinely 
Markovian process. 
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It is exactly at this stage where a possible justification of a linear response approach 
has to set in. In essence this requires two steps. First of all, if the 5t may be chosen to be 
sufficiently small on the collective time scale one may construct the K{Q, P,t;Qo, Po,to) 
by describing collective motion locally to harmonic order. Secondly, if the St is large enough 
on the microscopic scale the dynamics of the intrinsic degrees of freedom does not have to 
be followed in complete detail. Using such hypothesis it is possible to construct the form 
of K{Q, P, t] Qo, Pq, to) explicitly and to derive microscopic expressions for the individual 
transport coefficients. Moreover, this procedure even allows one to generalize Kramers' 
equation to include quantum effects, which show up, first of all, in generalized diffusion 
coefficients; for details see [2]. One step necessary in this direction is to interpret the 
K{Q, P, t; Qo, Pq, to) from above as a "propagator" for Wigner functions. 

It is clear that this locally harmonic approximation (LHA) is closely related to prop- 
erties which one expects to hold true for Fokker-Planck equations. Nevertheless, there are 
various ways to check that the goal set at the beginning is actually reached in the very 
end. For instance, it is possible to see a) whether the local propagators observe Marko- 
vian behavior, or b) whether or not the process is indeed "diffusive". As it turns out, the 
latter feature ceases to be given for unstable modes at low temperatures. We will not have 
time to go any further into these questions. We shall, however, be able to touch upon 
another condition for the LHA to be valid, the "smoothness" of the transport coefficients 
as function of the collective variables. 

2 Linear response theory for collective motion 

In the sequel let's suppose to be given a Hamiltonian H{xi,Pi,Q) for the nucleons' 
dynamics in a deformed mean field, with the deformation being parameterized by the shape 
variable Q, whose average (^H{xi,Pi,Q)) represents the total energy of the system E'tot 
(eventually including both the Strutinsky re-normalization as well as "heat"). The equation 
for average motion (EOM) for Q{t) can then be constructed from energy conservation. 
From Ehrenfest's equation it follows: 



All one needs to do to get the equation of motion for Q{t) is to express the average 






Following the scheme of the LHA one may expand 
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the H{Q) around any given Qq to have: 

1 /a^H 
H{Q{t)) = H{Qo) + m) - Qo)F + -m) - Qof ( QQ^iQo) ) (4) 

\ / Qo,To 

The effects of the coupling term {Q{t) — Qo)F may now be treated by linear response 
theory, exploiting as a powerful tool the causal response function x 

x{t-s) = Q{t-s)^tr (^MQo,To)[F\t),F\s)])=2iQ{t-s)x"{t-s) (5) 

Here, the time evolution in it) as well as the density operator are determined by 
H{Qq). The pqs is meant to represent thermal equilibrium at Qq with excitation being 
parameterized by temperature or by entropy. After some lengthy derivation one sees 
that the frequencies for local motion are given by the secular equation x(a;) + = 0, 
which actually determines the poles of the collective response Xcoii(<^) = x(<^)/(l + ^x(<^))- 
Different to common approaches but most important, in our case the coupling constant 
k appearing here is a derived quantity, given in the end by —k~^ = Id^H/dQ^) + 
X(0) - x^'^ = d^E{Q,So)/dQ^\^ + x(0), with x(0) being the static response, x""^ the 
adiabatic susceptibility and E{Q,So) the quasi-static energy at given Q and fixed entropy 
Sq. Finally, the transport coefficients for average motion can be introduced whenever it 
is possible to approximate the Xcoii('^) by an oscillator response function Xosc('^), in the 
sense of having 

iXcouiu;))-' S{F)^ ~ (Xosc((^))"' S{F)^ = {-Muj^-^ioj + C) 5{F)^ = -fe.ti^) (6) 

The /ext(<^) represents an "external" field which couples to our system through a term 
fex.t{t)F. Self-sustained motion corresponds to /ext(^) = 0, in which case the total energy 
must be conserved (according to (3)). As shown first in [5] (for the damped self-consistent 
case) the equation dE^Q^/dt = can be rewritten as 

-dt""^^'' = -dt [-^' + ) ^ ^^"^^^ = ""dt^ 

which correctly expresses the exchange between collective motion into heat. (The loi 
represents one of the possible (complex!) frequencies of the system, as determined from 
the secular equation). 
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3 Forced energy transfer to a system of independent particles 

As is clearly seen from (7), the friction force parameterizes that energy which is 
transferee! irreversibly to the intrinsic system. Let us study this feature within a simple 
model, with the simplifications consisting first of all in neglecting self-consistency. This 
means that we take a nucleus at given deformation Qq which is exposed to a time dependent 
external field at some fixed polarization. We may thus use a Hamiltonian of the type given 
in (4), where the F is chosen to represent this polarization, but where the last term on 
the right is neglected. The Q{t) — Qo = q{t) is then a truly external quantity, which shall 
be called q{t) in the sequel, and which is not subject to a subsidiary condition of the type 
k{F)t = q{t), which follows from (3) and (4). As another important simplification we 
will assume the H{Qo) to represent the ensemble of independent particles as given by the 
deformed shell model at zero temperature. 

Such a system has been studied in a series of papers which aimed at a new understand- 
ing of the physics of wall friction (see [3] and references given there) . The time dependence 
of the q{t) was assumed to be of the form q{t) = qo sin(f2t) and the system was followed 
for one period simulating the solutions of the Schrodiger equation numerically. 

Let us examine this problem within linear response theory. The energy transfered to 
the intrinsic degrees of freedom within one cycle may be evaluated from the following well 
known formula: 



The last expression is correct only for the truly periodic field. To derive it one needs to 
use the fact that the Fourier transform of the response function may conveniently be 
split into real and imaginary parts, x(c<;) = x'(^) + '>'X"{'^)j where (for real uj) x'i'^) is an 
even function and x"i^) an odd one. For this reason only the latter one survives after 
integrating twice over time. (For more details on these features see e.g. [2]). 

This result may be compared with those of [3], we simply need to identify AE^nt 
with their AE and calculate the total unperturbed energy Eq as the sum over single 
particle energies. However, one may as well go ahead and introduce already here a friction 
coefficient by the following reasoning. As the AEint measures the change of energy during 
one period of vibrations, one may simply divide by the length T = 2tt /Q. oi that period. 
In this way one gets: 





dEir^t _ T^qlx"{^) 



7(^)^2 



(9) 



dt {2tt/Q) 
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with V = q, qo = flqo and = {qo)'^/2. We have identified the friction coefficient as 

^(Q) ^ ^ ^ ^"(Q) (10) 

with the function being the so called relaxation function. Notice that the frequency- 

appearing here is the (real) one given by the external field. This is very different from the 
form indicated in (7). As mentioned there, the uji is the actual, complex frequency the 
collective motion has around the Qq. Incidentally, a form of the type x"{^)/^ may appear 
(for friction) even within the linear response formulation as described before, but only if 
the coupling between collective and intrinsic motion is treated perturbatively, see section 
3.3.2 of [2]. In case of small frequencies one may apply the so called zero frequency limit 
7(1^-0). 

In Figs.l and 2 we present numerical results for the quantity AE/Eq — AEi^t/Eo for 
the case of quadrupole excitations. They were calculated on the basis of our formula (9) 
but for the same system as in [3] , namely independent particles in a Woods-Saxon potential 
(of an un-physically large depth to decrease the escape probability). All parameters are 
chosen like there, which means that the 77 can approximately be written as r] ^ 0.02269^a;. 
As the most striking difference to the (quantal) results presented in Fig.l of [3], in our case 
we observe strong oscillations with r], which represent nothing else but the typical strength 
function behavior. (These functions are smooth in omega simply because we averaged the 
delta functions over an interval of 0.1 MeV). In both figures we show as the straight line 
marked with dots the result one gets in case that this energy transfer is calculated with 
wall friction. Apparently the latter result can be obtained at best after performing some 
averages. 

Indeed, it has been shown in [6] that the friction coefficient obtained within linear 
response theory (in the zero frequency limit) becomes close to the one of the wall formula 
after applying smoothing procedures in the sense of the Strutinsky method. This features 
goes along very nicely with the claim that wall friction represents the " macroscopic limit" , 
for a system of independent particles (for an extensive discussion of this topic see [2]). The 
same feature is seen here for the 7(ri) of (10). In Figs.l and 2 we present curves obtained 
from applying Strutinsky smoothing to the microscopic evaluations: For the dotted lines 
the averaging interval was 5 MeV, for the short dashed ones 10 MeV and for the long 
dashed ones 20 MeV. From Fig.l it is seen that and how smoothing leads to results similar 
to that of the wall formula. This calculation corresponds to the case where Qo stands for 
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a spherical deformation, the same situation which has been considered also in [6]. The 
one presented in Fig. 2 corresponds to a case where the unperturbed system has a sizable 
octupole deformation of a = 0.3. In this case the wall formula is not recovered (at least 
not the one corresponding to the spherical configuration used in the figure, whereas there 
could be some small dependence on deformation). We may say that similar results are 
obtained for vibrations of other multi-polarity. As an interesting feature we may note that 
"macroscopic" friction is smaller the more complex the microscopic strength distribution 
is. Wc would not like to speculate whether or not this fact is related to an increase of 
chaotic behavior of the nucleonic degrees of freedom. 

In [3] only forced vibrations around the sphere were considered. As can be seen from 
their Fig.l, the simulations of the Schrodinger equation indicate a straight behavior of 
the functional dependence of AE/Eq on rj somewhat below the wall formula. In such a 
non-perturbative calculation an average over a full cycle will differ in at least one respect 
from our procedure. Starting from formula (3) (which is a correct one), the (^F{xi,pi, Q)^ 
will no longer be a linear functional of q{s). If one still were expressing this quantity by an 
integral like ds x{i — s)q{s) the x{t — s) itself would have to be a complicated functional 
of q{s). Apparently, it behaves such that the average over the amplitude of oscillation (in 
the deformation degree of freedom) in the end leads to a linear relation with 77. It seems 
to us that this average may in a sense be considered analogous to our averaging in the 
spectrum, as is done in the Strutinksy method. As a matter of fact, experience tells one 
that averaging in energy e over an interval of = 10 MeV corresponds to an average in 
Q over a ■Jq^ — /ep ~ 1/4 — which corresponds nicely to the amplitude chosen in [3]! 

4 The influence of collisional damping on transport properties 

From the discussion of the last section it is clearly seen that for nuclear collective mo- 
tion it is not possible to justify a local friction force within the mere picture of independent 
particles. For such a model one has to employ averaging procedures of one kind or other. 
Moreover, we have observed that quite large intervals in the averaging parameters are in- 
volved if for the latter one chooses energy. This fact clearly hints to an inherent deficiency 
of the underlying model: At the excitations which are at stake in common experimental 
situations the picture of particles moving in a mean field without "collisions" does not 
applyl For this reason the notion of the H{Qo) to be simply given by the deformed shell 
model has been given up a long time ago whenever transport properties where calculated 
within the linear response approach (see [2] for a detailed discussion). Instead it was as- 
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sumed that the particles are dressed by self-energies having both real and imaginary parts: 
± ie, T) = T,'{(jJ, T) ^ f r(a;, T). The intrinsic response functions are then calculated 
after replacing the single particle strength Qk{oi>) = 2tt S{?ko — eu) by 

with the /U being the chemical potential. 

The I/Fq represents the strength of the "collisions", viz of the coupling to more 
complicated states. The cut-off parameter c allows one to account for the fact that the 
imaginary part of the self-energy does not increase indefinitely when the excitations get 
away from the Fermi surface. Both parameters are not known precisely, but from experi- 
ence with the optical potential and the efi^ective masses the following range of values can 
be given: O.OSMeF-^ < ^ < 0.06 MeV-^ and IhMeV <c< 30MeF. Neglecting the 
uj dependence of F and putting c — > oo these values lead to a average relaxation time for 
single particle motion Tint = h/T which is in accord with the estimate given in [7]. 

In Fig. 3 we present calculated along a fission path of ^^^Th for different temperatures 
(whose values are given in MeV). All curves except the one marked by triangles are identical 
to those of Fig. 13 of [8], where for the deformed shell model the Pashkevich code has been 
employed; details can be found in [8]. It is seen (i) that this ratio 7/M does not change 
very much with the collective variables as soon as T is of the order of 2 or larger, and (ii) 
that it increases with T (for reasons given below, the "heat pole" contribution has been 
removed in this calculation). For larger T the ratio is of the order as predicted by the wall 
formula (for 7) plus the one of irrotational flow for the inertia. The reason for this behavior 
is due to the fact that with increasing T the residual interactions become more and more 
important, with the two implications of a) smoothing out details of shell structure and b) 
making the microscopic mechanism of dissipation more effective. 

For T = 1 we have included a still preliminary result of an extension of our theory to 
the inclusion of pairing correlations. As expected the latter reduce the influence of shell 
effects, albeit details still will have to be clarifled further [9] 

5 The role of symmetries and the heat pole for nuclear friction 

Above it has been indicated that for the calculations presented in Fig.3 a particular 
contribution to friction was discarded. This shows up at flnite T and is related to an 
interesting quasi-static property which in turns is dominated by the influence of symmetries 
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(for a detailed discussion see [6, 2]). Let us demonstrate these features with the help of 
the zero frequency limit of friction. To sufficient accuracy the latter can be written as 

On the very right the correlation function has been introduced which is related to 
the response function by the famous fluctuation dissipation theorem (FDT) %x"{^) = 
tanh {?i(jj/2T) ip"{(jj). The definition of ip is similar to that given for x in (5), with two im- 
portant exceptions: The commutator is to be replaced by an anti- commutator ajid from the 
operator F one has to subtract its unperturbed average value (F) . The general microscopic 
expression for ip"{(jij) is 

= ^^0277(^(0;) + kV"(^) with V° = ^(x'^ - x(0)) (13) 

with the jii/j'^Lo) being regular at u> — 0. The is the isothermal susceptibility which 
measures how the (quasi-)static expectation value (F)'^^ changes with Q if the temperature 
is kept constant. The singularity at u; = is called "heat pole", in analogy to a similar 
pole in the density density strength distribution for infinite matter being responsible for 
heat diffusion there. A structure like that given in (13) is obtained only in the strict case 
of pure Hamiltonian dynamics (i.e. when the correlation function is formally calculated in 
the basis of exact eigen states) . Within our approximation of coUisional damping the heat 
pole changes like 



Both i^P as well as Ft have been calculated numerically in [8] . The result for Ft follows 
closely the following simple rule Ft ~ 2F(/x, T) 2T, which is valid over the very large 
range of temperatures of up to about 10 MeV. 

When applied to (12) one sees that the heat pole implies the following contribution 
to friction 

m •0° _ in 

^T 2r Ft 

Estimating x^ — x(0) simply in the independent particle model this component of friction 
turns into the one found first by Ayik and Norenberg within the model of DDD [10]. In 
[6] this form has been evaluated as function of temperature for all T. A slightly modified 
version is shown in Fig.4. The fully drawn line and the dashed one correspond the o7(0) of 



4/l lip 2% / rr. , ,\ 
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(15) with the c of (11) put equal to c = 20 MeV and c — > oo, respectively. The curve with 
the heavy dots corresponds to the contribution of the remaining part of the correlation 
function. As demonstrated in [6] (see also [2]), the distinction of the two contributions can 
simply be made in terms of the matrix elements of the (one-body) operator F with the 
shell model states. The o7(0) is solely to be associated to the diagonal elements. That 
they may lead to dissipation, nevertheless, (and thus to entropy production) is due to the 
effects of " collisions" . 

From Fig. 4 it is seen that at smaller T the o7(0) takes on very large values. They 
actually exceed several times that of wall friction (shown by the horizontal line), and they 
seem to be too large as required by experimental evidence (see [2]). The reason can be 
traced back to the following properties of static susceptibilities. Let us first rewrite the 
difference appearing in ■^^ as given by the left part of the following equation: 

^i^' = x^- x(o) = (x^ - x"-") - {x'" - x(o)) ^i^' = x^- X"'^ (16) 

The difference of the adiabatic to the isothermal susceptibility can be seen to be small in 
the nuclear case; it is proportional to the square of the cross derivative of the free energy 
with respect to Q and T, a quantity which for the system underlying Fig. 4 even vanishes. 
So the culprit must be the x^'^ — x(0)! However, this difference is known to vanish for 
truly ergodic systems, namely systems whose states are non-degenerate. As an additional 
condition one only needs to have a sufficiently narrow distribution of the energies of the 
excited states. 

Apparently, these conditions are not met for the case shown in Fig. 4, and most likely 
both of them are violated. First of all, the microscopic evaluation of the matrix elements 
is dominated by the model of independent particles, with all the many degeneracies ap- 
pearing there. Secondly, applying the canonical distribution to parameterize the thermal 
excitations of a nuclear system the spread in energy is exaggerated artificially. On the 
other hand, there is little doubt that the true compound configurations will remove these 
spurious contributions. First of all, because a consideration of the compound states will re- 
quire a more correct treatment of the residual interactions. In this way the many artificial 
degeneracies of the deformed shell model will be removed. Secondly, thermal excitations 
will have to be treated on the basis of the micro canonical ensemble. To simply simulate 
these effects one may just apply the reduction indicated in (16) to remove the unphysical 
contribution from the heat pole. 



10 



6 Dissipation within Landau theory 

As seen above, the friction coefficient tends to decrease with T at larger temperature. 
This feature is evident for the component o7(0) (see (15) and Fig. 4), but as discussed in 
[6] it wiU eventually hold true also for the other component (see also [11]) under certain 
circumstances (like approximating the imaginary part of the self-energy in "common" 
relaxation time approximation (with c = oo)). 

Such a behavior with T reminds one of the two body viscosity of hydrodynamics. 
In [12] a model has been suggested in which the intrinsic dynamics is described by the 
Landau- Vlasov equation and where the finiteness of the system exhibiting shape dynamics 
is introduced through special boundary conditions. In Fig. 5 we present a calculation of 
the friction coefficient (as function of T) for quadrupole vibrations about a sphere, done 
within an extension of this model. The dashed and the fully drawn lines correspond to the 
hydrodynamical limit, for two different choices of the parameter c entering the relaxation 
time used in the collision term of the Landau- Vlasov equation. The squares correspond 
to contributions from different peaks in the correlation function, where the full ones are 
suposcd to correspond to the analog of the "heat pole". Similarities with the behavior 
shown in Fig. 4 are evident. So far however, it is yet unclear exactly where the contribution 
comes from, solely from the difference ~ X^'^ ^is it should, or whether also in this model 
there is the spurious effect coming from a non- vanishing difference x^'^ ~ x(0)- Further 
studies are under way. 

7 Summary 

To describe collective motion as a Markovian transport process one needs to be able 
to define transport coefficients which vary smoothly with the macroscopic variables, which 
by the way have to include the parameter which measures the thermal excitation. We 
have demonstrated that such a condition can hardly be fulfilled within the picture of the 
deformed shell model. On the other hand, we have shown that residual interactions may 
do the job, the better the larger the thermal excitation. At present the situation is less 
clear at smaller temperatures. Whether or not pairing alone will do is currently under 
investigation. It may well be, however, that even in this regime one may want to include 
more of the configurations as given by the nuclear compound model. 
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Figure captions 

Fig.l: Average energy transfer to a spherical system of independent particles by an external 
quadrupole perturbation, calculated within linear response. Otherwise the same picture is 
adopted as in [3]. 

Fig. 2: Same as for Fig. 2 but for a system with octupole deformation. 

Fig. 3: Ratio of friction to inertia along the fission path of ^^^Th (for details see text). 

Fig. 4: The contribution of the "heat pole" to friction for collective quadrupole oscillations 
of a system of particles in a square well potential (see text). 

Fig. 5: Friction for quadrupole oscillations calculated from a Landau- Vlasov approach to a 
finite nucleus. 
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